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1. Introduction 

The gauge-gravity duality conjecture, which is a generalization of the AdS/CFT 
correspondence, relates theories of quantum gravity living on curved spacetimes 
and quantum gauge field theories without gravity defined on the boundaries of 
such spacetimes. This conjecture provides very promising directions for the study 
of quantum gravity and black holes. It is possible to describe certain black holes 
in terms of the worldvolume theories of the D-branes that compose them. These 
theories are the maximally supersymmetric Yang-Mills (SYM) theories in various 
dimensions, taken in the’t Hooft limit and at finite temperature. They are strongly 
coupled in the regime in which they describe string theory black holes. Solving 
these thermal gauge theories at strong coupling would allow one to directly study 
the quantum properties of the dual black holes, including their thermodynamic 
aspects. 

The AdS/CFT correspondence relates maximally supersymmetric Yang-Mills 
theories with gauge group SU{N), taken in the large N ’t Hooft limit to certain 
closed superstring theories in the near horizon region of N coincident D-branes. The 
simplest case, which is the main focus of this review, is the one in which the gauge 
theory is the SU(N) supersymmetric Yang-Mills quantum mechanics and its dual 
theory is type IIA string theory describing N DO-branes. 

In this review, we focus on the recent developments in the gauge-gravity quality 
conjecture emphasizing the results from studying the strongly coupled supersym¬ 
metric gauge theories in one and two dimensions. To study these gauge theories 
at strong coupling directly we need to formulate them in a non-perturbative set¬ 
ting, and this can be achieved by regularizing them, in a consistent manner, on a 
spacetime lattice. We also discuss the lattice regularizations - naive and manifestly 
supersymmetric - of maximally supersymmetric Yang-Mills theories and the sim¬ 
ulation algorithms needed to compute the interesting observables of the thermal 
gauge theories at strong coupling. 

This review is organized as follows. In section [2] we briefly discuss the quantum 
gravity corrections to the internal energy of a black hole from low energy string 
theory calculations. We write down the internal energy of the quantum black hole 
as a function of its temperature with the anticipation that we can reproduce such 
a behavior from non-perturbative simulations of the dual thermal gauge theory. In 
sec. [3] we write down the action of the sixteen supercharge Yang-Mills quantum 
mechanics theory - the BFSS matrix model - that takes part in the gauge-gravity 
duality. The BFSS model serves as the simplest setting in which we can probe the 
duality conjecture and related quantum gravity corrections. In sec. 13.II we discuss a 
twisted version of the BFSS matrix model which shows compatibility with the lattice 
regularization by preserving a subset of the supersymmetry charges at finite lattice 
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spacing. In sec. |4]we write down the lattice actions for maximally supersymmetric 
Yang-Mills theories that would allow us to explore the strong coupling regimes 
of these theories. In sec. 14.11 we write down a naive lattice action for the sixteen 
supercharge gauged quantum mechanics while in sec. l4.2l we write down a manifestly 
supersymmetric lattice action of the same theory. In sec. [5] we provide details of the 
simulation algorithms used to probe the strongly coupled regimes of these theories. 
We provide results from the recent lattice investigations in sec. [SI which support the 
predictions on the quantum gravity black holes from the gravity side of the duality. 
In sec. 13 we discuss a one-parameter deformation of the BESS matrix model, known 
as the BMN model. This model is particularly interesting from the simulation point 
of view as it does not suffer from the thermal divergence of the partition function 
at smaller number of colors. In sec.|8|we provide a brief outlook on the supergravity 
system with Dl-branes and associated black hole - black string phase transitions. 
In sec. we discuss supersymmetric lattices for the two-dimensional maximally 
supersymmetric Yang-Mills. In sec. l8.2l we provide results on the Gregory-Laflamme 
phase transitions from lattice simulations of the strongly coupled gauge theories. 
We conclude in sec. |9|with the current status on validating the duality conjecture 
from the gauge theory side of the correspondence and possible future directions. 


2. DO-branes and Quantum Gravity Corrections 

String theory contains solitonic objects known as D-branes. In general, in a full 
string theory, there are modes that propagate in the bulk and also modes that 
propagate on the solitons. It is possible to define a limit of the full theory in which 
the bulk modes decouple from the modes that live on the solitons. This is typically 
the low energy limit of the theory, and in this limit, the D-brane theory becomes 
a supersymmetric Yang-Mills (SYM) theory. Since D-branes can carry mass and 
charge we can find supergravity solutions carrying the same mass and charge. Su¬ 
pergravity solutions are valid as long as the curvatures are locally small compared 
to the string scale Ig. When we have a system with large number N of D-branes the 
curvature is small and we can trust the supergravity solution even at sub-stringy 
distances. Supergravity solutions corresponding to (O-l-l)-dimensional SYM theory, 
which we are interested in this review, are black DO-brane solutions. 

We study DO-branes in the field theory limilPt^ 

9ym = “ fixed, a =11 ^ 0, (1) 

47r"^ a'2 

where Qs is the string coupling and gyM the Yang-Mills coupling. We keep the 
energies fixed when we take the above limit. In our case it implies that the theory 
decouples from the bulk since the ten-dimensional Newton constant goes to zero. It 
also suppresses the higher order a' correction^ in the action of the theory. 


'They come from the effects due to the strings having finite length. 
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When we take the limit given in Eq. (HD we are interested in the finite energy 
configurations in the field theory. This corresponds to finite Higgs expectation values 
for some helds of the theory. We are therefore considering the limit 

U = — = fixed, a' —)> 0, (2) 

a 

with r representing some position at which the DO-brane is sitting, and it gives a 
Higgs expectation value U to some fields of the theory. 

We will also consider near extremal configurations which correspond to the de¬ 
coupled field theories at finite temperature. On the supergravity side we start from 
a near extremal black DO-brane solution and we take the limit Eq. ([T]) keeping the 
energy density of the brane finite. In this limit the metric near the horizon has the 
following forrrPl^ 

(~ 27rXo^^A ^ 27rV240^^A {u-i^ + ^ , (3) 

with the function 

/([/) = I-(4) 

and representing the line element of an eight-dimensional unit sphere. The 
event horizon of the black hole is located at U = Uq. The parameters A and Uo are 
related to the mass and charge of the supergravity black hole. 

The black hole described above appears as a quantum system with temperature 
T for a distant observer. It undergoes thermal evaporation through black body 
radiation. The Hawking temperature characterizing this is obtained by a standard 
calculation 


As 167r3-^157r VAs / 


( 5 ) 


The Bekenstein-Hawking entropy of the black hole is evaluated by the area A of 
the horizon in the Einstein frame a^ 


S 1 A y26ig2^14N 5 / y N I 


The gauge-gravity duality tells us that T and S should be identified with the 
temperature and entropy of the dual gauge field theory. We can make use of the 
first law of thermodynamics, dE — TdS, to obtain the internal energy E = E{T) 
of the black hole. Normalizing by As we have 


E 


9 

14 




= 7.41 


( 7 ) 


where E and T are now dimensionless variables. 

It is possible to add quantum corrections to Eq. ©■ The back reaction of the 
Hawking radiation can contribute as non-local quantum corrections to the super¬ 
gravity action. The magnitude of such corrections can be estimated by assuming 
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the Stefan-Boltzmann law and using the black hole geometry given in Eq. ([3]) . The 
calculations show that the energy loss due to Hawking radiation is suppressed by a 
factor of 0{ll^) and it indeed vanishes in the low energy limit, Ig —>■ 0, considered 
here. 

One could also add another correction to Eq. 0 due to the quantum gravity 
effects arising from short distances. One can obtain such a correction by calculat¬ 
ing perturbatively the scattering amplitude involving four gravitons as asymptotic 
states in superstring theory. This amplitude is coming from an intermediate process 
of strings being pair-created and pair-annihilated near the horizon. Such quantum 
corrections can be implemented into the supergravity action by adding quadratic 
terms of the Riemann tensor at the leading order as shown by Gross and Witten.l^ 
Including such quantum corrections one obtains the metric near the horizon as 


where 




ds^ = 


240^5^ /C7j ^ 


( 8 ) 

(9) 


and hi, i = 1,2 and /i are functions of U/Uq. The functions Hi, 2 {U) and Fi{U) 
can be determined uniquely to the order I/V2 by imposing appropriate boundary 
conditions as shown by Hyakutake.^^ The metric given in Eq. ([8|) reduces to the one 
in Eq. ([3]) when iV —>■ 00 , which corresponds to the limit of classical gravity. It is 
seen from Eq. that the position of the event horizon is slightly shifted inward 
due to quantum corrections as 


u = ( 11 ) 

where fi(Uo) is some positive constant. 

In the modified background geometry given in Eq. (|S]) one can proceed to calcu¬ 
late the Hawking temperature in a standard manner. We need to use Wald’s entropy 
formula, which is applicable to the present case with higher curvature terms, in order 
to evaluate the entropy of the quantum black hole. The temperature and entropy of 
the black hole can be identified with the corresponding ones in the dual gauge the¬ 
ory with the contributions coming from the 1/N'^ corrections. The internal energy 
E in this case can also be obtained directly from the first law of thermodynamics. 
Thus we arrive at 


E 


= 7.41 - 5.77 T 


0.4 


iV2’ 


( 12 ) 


where the second term represents the quantum gravity corrections at the leading 
order. 
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We have ignored the a' corrections in the above analysis. One can include these 
corrections to Eq. (nm as has been done in Ref. [5] at V = oo. This gives the 
following result for the internal energy of a quantum black hole 

^iSgravity = (7.41 + A + . . . ) 

+ (-5.77 rO-4 + B + ...)^ + o , (13) 

where the constants A and B are not yet known. We have denoted the internal en¬ 
ergy by Egravity in the above expression since this result is coming from calculations 
on the gravity side of the duality. We will denote the internal energy obtained from 
the gauge theory side of the correspondence as Egauge- The power of T for each 
term in the above result can be determined from dimensional analysis using some 
known results in superstring theory.^^ 

3. The BFSS Matrix Model 

The BFSS matrix modeP^ is the maximally supersymmetric matrix quantum me¬ 
chanics, which is obtained by dimensionally reducing four-dimensional maximally 
supersymmetric Yang-Mills theory to one dimension. This model is conjectured to 
describe a discrete light-cone quantization of M-theorjff^l on a circle. It also de¬ 
scribes, with suitably re-identified parameters, the low energy dynamics of a col¬ 
lection of N DO-branes of type IIA superstring theory.l^^ This theory also plays 
the role as a matrix regularization of the light-cone action for the 11 -dimensional 
super-membrane. (See [13] for a review.) 

The theory contains a U{N) gauge field Aq, nine scalars Xi {i = 1, ■ ■ ■ ,9), which 
are N x N hermitian matrices, and 16 fermions ip a (a = 1, • • • , 16). The continuum 
Euclidean action of the theory has the form 

^ -h ^IpaDrlpa 

+ ^'0a(7i)a/3[^i,V'/3]}, (14) 

where A = NgyM is the ’t Hooft coupling and 7 ^ are real symmetric matrices 
satisfying the nine-dimensional Clifford algebra, { 7 i, 7 j} = 

The covariant derivative is defined as 

Dt' = dr ■ +i[AQ, ■ ]. (15) 

The theory is invariant under the following set of supersymmetry transforma¬ 
tions 

^^Xi — (7^)q:/3'0/3 ; 


(16) 
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where 



(17) 


and are sixteen fermionic parameters. The invariance of the action Eq. (I14p under 
the set of transformations Eq. (1161) can be shown using the Fierz identity for the 
fermions and the Jacobi identity for the scalars. 

The supersymmetry charges Qa of the theory are the generators of the infinites¬ 
imal transformations, 


(18) 


— ^aQa ■ 


In Sec. 63D we study a lattice version of this theory, which is obtained through 
a naive lattice regularization. In Sec. (14.21) we look at another lattice discretization 
of this theory, which is manifestly supersymmetric; it preserves two of the sixteen 
supersymmetry charges on the lattice. 

3.1. Topological Twist of the BFSS Matrix Model 

The BFSS matrix model can be expressed in a convenient form, which is compati¬ 
ble with lattice regularization, after relabeling its fields and supersymmetries. The 
lattice theory obtained through such a prescription is manifestly gauge-invariant 
and preserves two or eight supersymmetry charges exact, depending on the reg¬ 
ularization scheme, at finite lattice spacing. This method of relabeling the helds 
and supercharges is known as topological twisting. The twisting process is done on 
its four-dimensional parent theory, the A/” = 4 SYM theory. The four-dimensional 
theory has 5'0(4) Euclidean spacetime rotation group and SO{6) internal rotation 
(R-symmetry) group. We introduce a new rotation group 5'0(4)' by taking the 
diagonal subgroup 


50(4)' = diag subgroup (50(4)Lorentz x SO{A)r) . 


(19) 


The theory obtained this way is called twisted theory. The fields and supersym¬ 
metries of the original theory are then expressed in terms of the twisted Lorentz 
rotation group. 

We dimensionally reduce the four-dimensional twisted theory to one dimension 
to obtain a twisted version of the BESS matrix model. This model can be easily 
transported on to the lattice by preserving gauge symmetry and a part of the 
supersymmetry. In this section we discuss the twisted theory that leads to a lattice 
action proposed by Sugino. In Ref. [Ml [151 E] Kadoh and Kamata performed a 
series of lattice simulations of this theory and examined the gauge-gravity duality 
in this system. We note that there also exist a complementary lattice discretization 
of the 16 supercharge theory, based on the idea of orbifolding, first proposed by 
Kaplan and Unsal.l^lJ^ This construction preserves 8 supercharges at finite lattice 
spacing. It would be interesting to carry out a detailed nonperturbative numerical 
study in the system prescribed by such a regularization scheme. 
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Let us briefly discuss the twisted version of the BFSS matrix model. The trans¬ 
formation laws of the fields and supercharges under the twisted Lorentz group are 
different from those of the untwisted theory. We can construct two linear combina¬ 
tions of supercharges from the original sixteen supercharges, which we label as Q±. 
These two supercharges, which we call as twisted supercharges, are scalar fermions 
and nilpotent up to a gauge rotation. After twisting, we can express the action of 
the theory in a topological (3±-exact form. Then the (5±-invariance of the action 
follows from the nilpotency of Q±. 

We change the field variables of the original theory to those of the twisted theory 
in order to define the Qi-exact action directly from Eq. (HI. The bosonic twisted 
fields are denoted by A^,Bi,C and (/)±, with = 0,1,2,3 and i = 0,1,2, and 
the the fermionic twisted fields are denoted by r]±,ip±fj, and x±i- Under the four¬ 
dimensional twisted Lorentz transformations the fields and ip±^ transform as 
vectors, the fields C, 4'± and rj± remain unchanged (they transform as scalars), and 
the other bosons and fermions of the theory transform as self-dual tensors. The 
fields Bi and x±i three independent components of the bosonic and fermionic 
tensors, respectively. In one dimension we associate the indices /x, i to the names of 
the fields. 

The two nilpotent twisted supercharges Q± can be constructed as linear combi¬ 
nations of the supersymmetry charges of the original theory in the following way 


Q+ = "^(*55 + fQia), 
Q- = + fQg), 


( 20 ) 

( 21 ) 


where Qi denotes the i-th component of the sixteen-component supersymmetry 
charge. 

The twisted fields in the bosonic sector are related to the untwisted fields the 
following way 


a^=1,2,3 (22) 

V 4 = -B 2 V 5 = (23) 

Ag = -So ^7 = \C (24) 

^8 = ^(<(>+- <('-) Xg = + (j)_) (25) 

with Aq unchanged. The fields A^, Bi and C are hermitian and ((/'+)^ = —</'-■ 


The fermionic twisted fields rj±,ip±^ and x±i ^re related to the components of 
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the untwisted sixteen-component fermion field 'ifja the following way 


(i’l \ 


/ V'-o + |?7+ \ 

V'2 


i-i - iX+2 

V'3 


i -2 + ix+i 

V'4 


i-3 - *X+0 

V'S 


i+0 + ^v- 

V'e 


i+i + *x-2 

V'7 


i +2 - ix-i 

V'S 

1 

i+3 + *X-0 

V'g 

“71 

-i{i-o - §»7+) 

i’lo 


-i{i-i +ix+ 2 ) 

i’ll 


-i{i -2 - ix+i) 

ii2 


-i{i-3 + fx+o) 

ii3 


-i{i+o - §77-) 

in 


-i{i+i - ix- 2 ) 

iib 


-i{i+2+iX-i) 

\iw/ 


\-i{i+3 - ix-o)/ 


We use appropriate set of gamma matrices (see Ref. [1^]) to connect the two sets 
(twisted and untwisted) of the field components. 

It is straightforward to obtain the (5±-transformation laws on the twisted fields. 


They are given by 

Q±Afj, = (27) 

Q±ip±ij. = -iDf,(j)± (28) 

Q±Bi = x±i (29) 

Q±X±i = [Bi,4‘±] (30) 

Q±c = r,± (31) 

Q±V± = [C,cl^±] (32) 

Q±4’t ~ (33) 

Q±Vt ~ (34) 

Q±^± = 0 (35) 

Q±XTi = l[C,B,]±H, (36) 

= (37) 

Q±Hi = ± ^[xzpi,(?!>±] + ]f[X±i,C] + ]fB^,'q±^ (38) 

= ± + \[x±p..C] - (39) 

where Hq is the covariant derivative and with ^ = 1 , 2 ,3, are the commutators 
Dfj_- = ■ ]. We have used the four index notation for (/r = 0,1,2,3) 
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since the twisted theory is originally defined in four dimensions. We have also in¬ 
troduced seven auxiliary fields Hi and to make the Q± transformations satisfy 
the following relations 

Ql=td^^, (40) 

{Q+,Q-} =-iSc, (41) 

where the right-hand sides correspond to gauge rotations with field dependent gauge 
parameters. 

After relabeling the variables of the original theory we write down the action of 
the twisted theory in a Q±-eiiSLCt form 

N 

S = ^Q+Q- A, (42) 

where the “gauge fermion” A is 

A = y dr Tr I - 2iB^(^F^3 + ^eijkFjk^ + Bk] 

-V’-TmV'-m + X+iX-^ - , (43) 

with Fqi = DoAi and Fij = i[Ai,Aj\ with i,j = 1,2,3. 

We see that the twisted action is trivially invariant under the Q±- 
transformations, without invoking the Leibniz rule, owing to the nilpotency of the 
twisted supersymmetry charges. The other fourteen supersymmetry charges and the 
global S'0(9) symmetry are also exact symmetries of the continuum twisted action 
since we have only performed a change of variables of the original theory to obtain 
the twisted theory. 

4. Lattice Actions 

There are two ways to discretize the sixteen supercharge quantum mechanics model 
on the lattice. One is a naive discretization of the continuum theory and the other 
is a manifestly supersymmetric lattice action. The naive lattice action follows the 
standard rules of discretization of a gauge theory. The resulting fermion determinant 
can be complex in general at finite lattice spacing, which can make the simulation 
problematic. While the naive lattice action preserves no supersymmetry, we can 
argue that since the quantum mechanics is free from UV divergences given that the 
gauge and global symmetries are preserved by a UV regulator, all the supersym¬ 
metries will be restored in the continuum limit. The twisted version of the theory 
is manifestly supersymmetric and gauge-invariant on the lattice. It preserves two 
supersymmetry charges exact at finite lattice spacing while the other fourteen su¬ 
percharges are broken on the lattice. Supersymmetry is broken softly in both the 
lattice constructions by thermal boundary conditions and also by the addition of 
mass terms for the scalars to stabilize the simulations. 

We discretize the theory on a Euclidean time circle consisting of M number of 
sites, n = {0, • • • , M — I}. Thermal boundary conditions correspond to taking the 
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fermions anti-periodic on the Euclidean time circle. The temperature of the system 
is T = 1//3, with /3 = aM and a denotes the lattice spacing. 


4.1. Naive Lattice Action 


The continuum Euclidean path integral of the sixteen supercharge theory takes the 
following from after the fermions are integrated out 


Z = J DADX Pf(7W)e"^^, 


(44) 


where Pf(A^) denotes the Pfaffian of the fermion operator A4. The bosonic action 
is 




dr 




and the fermion operator 


M = 'JrDr - ■ ]. 


(45) 


(46) 


We take 7 ,-, 7 ^ in the Euclidean representation of the Lorentzian Majorana-Weyl 
gamma matrices, obeying { 7 i, 7 j} = 2dy for index v = {r, i}. The Pfaffian is in 
general complex,^ giving rise to a potential sign problem. It is thus important, in 
principle, to include the phase of the Pfaffian in the Monte Carlo simulation. If the 
phase of the Pfaffian is very nearly real and positive we can ignore the phase and 
proceed with phase quenched simulation^. 

We also need to add mass terms to the lattice action in order to regulate the 
flat directions associated with the scalars in the theory. The classical potential, 
^ Tr {[Xi, Xj]'^), has flat directions, that is, any set of matrices which are mutually 
commuting have zero energy. These flat directions can be removed by adding mass 
terms of the form 

r/5 ^ 

( 47 ) 

i=l 

with /r a mass parameter. In general one simulates the theory at various values of 
the mass parameter and then extrapolates to vanishing mass limit. 

We discretize the continuum model as. 


Af ^ 1 

^M = ^f dry]Tri/i2(v2), 




M-1 


n—0 


and 


XI mn — 


0 

-(D.) 


0 




(48) 


(49) 


^See Refs. |20l 1211 \22\ for a recent work on the sign problem in 4 and 16 supercharge Yang-Mills 
theories in two dimensions. 
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where we have rescaled the fields Xi^n and ipa,n by powers of the lattice spacing to 
render them dimensionless. The dimensionless lattice mass regulator fj,Q is related 
to the continuum mass as 

Mo = /3 m> (50) 

and it is this dimensionless fig that parameterizes how much of the weakly-coupled 
divergent region is allowed by the mass regulator.!^ 

We define the lattice covariant difference operator the following way 

- UnW,,n-lUl (51) 

and we have introduced a Wilson gauge link field [/„. Note that the fermionic 
operator M is free of doublers and manifestly antisymmetric on the lattice. We 
have chosen a twisted Euclidean representation for the gamma matrices in order to 
obtain this antisymmetric fermion operator after taking finite differences. We have 
taken a representation of gamma matrices where 

In a super-renormalizable theory only a finite number of divergences can arise 
in the lattice regularized theory. Such divergences may induce a RG flow away 
from the maximally supersymmetric continuum limit of our interest. However, one 
would naively expect that for quantum mechanics there will be no divergences. Thus 
we expect that the lattice action is finite and hence will flow without fine-tuning 
to the correct maximally supersymmetric continuum theory as the lattice spacing 
is decreased. In Ref. [53] this was tested by simulating the theory with periodic 
boundary conditions and testing whether the path integral, the Witten index, was 
dependent of /3 in the continuum limit. 

4.2. Manifestly Supersymmetric Lattice Action 

It is straightforward to discretize the topologically twisted theory given in Eq. (I321). 
We place the adjoint scalars and fermions on the sites and the lattice gauge field 
Un S SU(N) on the links. The covariant difference operator on the lattice acts the 
following way on a generic field ipn 

— Uji^n+lUn (5fi) 

and it transforms as (pn itself. 

By replacing the integral and the covariant derivative in the twisted action Eq. 
(H5]) with the summation over the lattice sites and the covariant difference operator, 
respectively, we obtain the following lattice action (We call this the Sugino lattice 
actiorP^ as it is the dimensional reduction of the four-dimensional lattice theory 
constructed by Sugino.) 

N 

S = Q+Q-— A, 


( 54 ) 
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where 


M-1 

A = XI ( 

n—0 

“t“ ^ ^k\ 

+X+*X-* - 

(55) 

with 



Foi = D^Ai and 

Fij = ^[^2, h j ~ 

(56) 

and we have denoted 

Aq = 

(57) 


as the dimensionless ’t Hooft coupling. The continuum limit of this theory is ob¬ 
tained by taking the limit Aq —^ 0 with A fixed. 

The Q± transformations on the twisted lattice fields take the following form 


Q±[/ = 

(58) 

Q^Afj^ = 

(59) 

Q±V’±/x = -iDfj,(j)± + iSfj,oilj±oilj±o 

(60) 


(61) 

Q±x±i = [BiA±] 

(62) 

Q±C = ?7± 

(63) 

Q±-n± = [C, (j)±] 

(64) 

Q±'t’T ~ Vt 

(65) 

Q±Vt = 

(66) 

Q±(j)± = 0 

(67) 

Q±XTi = ± Hi 

(68) 

Q±‘0=F/j. — zb o} 

(69) 

Q±Hi = ± + ^[x±i,C] + 

(70) 

Q±Hf, = ± + ^[x±y,C] - ^Df,r]±]\ 



(^^bP±o,DoC ±2iH {tp+o,tp-o}] + (71) 

They are similar to the corresponding continuum transformations with the exception 
of the fields with index ft = 0. The terms containing S^o stre higher order corrections. 

The supersymmetry charges Q± satisfy the following set of relations on the 
lattic^^ 


Q± — , 

{Q+,Q-} = i6c, 


(72) 

(73) 






October 12, 2015 0:50 WSPC/INSTRUCTION FILE REVIEW 


14 Anosh Joseph 

which resembles the corresponding set of relations in the continuum. Such a resem¬ 
blance does not come as a surprise if we recall the BRST transformations. In the 
BRST transformations, once the transformation rule for a lower dimensional field 
is given, then that of the higher dimensional field is uniquely fixed from the nilpo- 
tency property. The same can be applied on the lattice, once the Q^-transformation 
laws of the lower dimensional fields, Q±U,Q±Bi, etc., are given, then the others, 
Q±'ip, Q±X±j s.tc., are uniquely determined from the above nilpotency equations. 

The lattice action has -invariance and gauge-invariance. We note that the 
remaining supercharges and the SO(9) global rotational symmetry are not pre¬ 
served on the lattice. The lattice action has, at least, an SU{2) global symmetry 
corresponding to the interchange of (5+ and Q-. The lattice Q± transformations 
reproduce the corresponding set of continuum ones and the lattice action reproduces 
the correct continuum action in the naive continuum limit. 

After performing (5±-variations on the lattice action Eq. (l5^ we find that it has 
a four-Fermi term. 



(74) 


which is of the order of the cut-off. The four-Fermi interaction term is not suitable 
for numerical simulations. One could express the four-Fermi term as a fermion 
bilinear form, as in the NJL model, by introducing an auxiliary scalar field a. The 
four-Fermi term becomes 



(75) 


n 


It reproduces the original four-Fermi piece once a is integrated out. 

5. Lattice Simulations of Sixteen Supercharge Yang-Mills 

The maximally supersymmetric Yang-Mills theory can be simulated on the lattice 
using rational hybrid Monte Carlo (RHMC) algorithm introduced in Ref. [26]. We 
are interested in numerically computing operator expectation values {O) through a 
discretized version of the path integral 



(76) 


where the partition function is 



(77) 


and $ and denote, respectively the set of bosons and fermions in the theory. 
The discretized path integral is evaluated through the importance sampling Monte 
Carlo algorithm, which we briefly discuss below. 
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We first integrate out the fermions in the theory and this leads to a Pfaffian. 
The partition function becomes 

Z = j D<^ Pf(7W($)) (78) 

The Pfaffian of the theory can be complex in general, 

Pf(7W) = |Pf(7W)|e*“, (79) 


with a phase a and it can affect the stability of numerical simulations. Investigations 
on the sixteen supercharge lattice theory show that this model does not suffer from 
a severe sign problem - the Pfaffian is very nearly real and positive.l^^®^ Thus we 
can ignore the phase of the Pfaffian and use the absolute value 

|Pf(7W)| = |det(7W)p = det (TW^TW)^ . (80) 


The next step is to introduce a set of bosonic pseudo-fermions ^ and express the 
absolute value of the Pfaffian 

|Pf(7W)| = y’z?^^£>Cexp(-C^(7Wl'7W)-ie) , (81) 

where ^ is a I6(V^ — l)-component vector (for SU{N) gauge group) on each lattice 
site. 

At each site n the pseudo-fermion field is set to 

= (82) 


with r]n being 16(V^ — I)-component random vectors drawn from a Gaussian dis¬ 
tribution. 

In order to simulate the theory we need to approximate the rational exponent 
of the squared fermion operator with a partial sum 

(7Wl7W)"i = Qfo+ TTTTirTTT’ (^3) 


A4^Ai+Pi 


where the degree P, and the coefficients ai, Pi known as amplitudes and shifts, 
respectively depend on the range of the eigenvalues of Ai^Ai and the accuracy 
of the approximation. We can compute the optimal shifts and amplitudes offline, 
through the Remez algorithm,!^ by minimizing the relative approximation error 
within a given spectral range [Aiow, Ahigh], with Aiow > 0 of the squared fermion 
operator. We need to add more terms to the partial fraction approximation if we 
want to keep the errors negligible across larger spectral ranges and this in turn 
results in increased computational costs. We must demand, at a minimum, that the 
smallest and largest eigenvalues of AVAi lie within the spectral range we use, that 
is. Alow < Amin ^ Amax < Ahigh- We should also monitor the extremal eigenvalues 
of the squared fermion operator, like the Pfaffian, during the RHMC evolution. 

Following the standard procedure, we then introduce fictitious momenta 
drawn from a Gaussian distribution, that are conjugate to the bosonic coordinates 
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<I> and evolve the coupled system using molecular dynamics (MD) evolution along 
a trajectory of length r in a fictitious simulation time t. The effective Hamiltonian 
has the form 

i? = ^ E Tr (4J + SB{^n) + (84) 

n 

We keep the pseudo-fermion field fixed along the trajectory as it already has a 
proper distribution. Since the MD evolution involves inexact integration of Hamil¬ 
ton’s equations, which does not conserve the effective Hamiltonian, we need to 
impose a criterion to accept or reject the new set of field configurations at the end 
of the trajectory. The new field configurations produced by the MD evolution are 
accepted or rejected using the Metropolis-Rosenbluth-Teller test, with a probability 

p = min(l, (85) 


where AiL is the change in effective Hamiltonian after one MD evolution. The 
original configuration values are restored back if the new configurations are not 
accepted and the evolution is repeated from a fresh set of Gaussian distributions 
for the momenta. The test makes the algorithm exact by stochastically correcting 
for integration errors as long as the MD integration scheme we use is area preserving 
(symplectic) and reversible (symmetric). 

Each step of the MD evolution requires solving the following set of equations 


It 

St 


SH 

Sp<!> 


m _ ssb 

(5<i) (5$ 




( 86 ) 

(87) 


The gradients of the fields with respect to the pseudo-fermion action (the 
fermionic forces) have the form 


_S_ 

(54) 




2=1 


54 ) 


X 




( 88 ) 


The quantities + Pi) can be efficiently determined by a multi-shift 

conjugate gradient (CG) inverter.!^ 

The Pfafhan phase may be re-incorporated in the expectation value of an ob¬ 
servable O by re-weighting as 


{O) 


E^(e*“) 


(89) 


where the sum runs over all members of the phase quenched ensemble. 
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6. Quantum Gravity Results from the Lattice 

In this section we discuss the results from Monte Carlo simulations of the maxi¬ 
mally supersymmetric Yang-Mills theory at finite temperature. The most interesting 
quantity to compute is the internal energy Eg^uge of the system. It has the form 

^L'gauge = ~ ~ Sq) , (90) 

where Sq denotes the extensive zero-point energy contribution and it is given by 

So = ^M{N^-l). (91) 

The other observables of interest are the expectation value of the Polyakov loop 

(P) = ^Tr [/, (92) 

with the holonomy matrix U dehned as 

U = Pexp \ i (f (irAo(T) I , (93) 


with P denoting a path-ordered product and the “extent of space” 

(R^) =f dr Tt {X,)^y (94) 

The expectation value of the Polyakov loop plays the role of an order parameter 
for the confining-deconfining phase transition in the thermal theory. It is non-zero 
at all temperatures in the supersymmetric theory. This implies that there is no 
phase transition in the supersymmetric theory as predicted by gauge-gravity corre¬ 
spondence.!^^ 

One could also check if the Polyakov loop can be fitted to the form 

(P)=exp(-|; + 6), (95) 

which is a characteristic behavior in a deconfined theory. 

The instability related to the flat directions associated with the scalar eigenvalues 
can be probed by studying the divergence of (R^)- In Ref. [32] it has been shown, 
for numerical simulations based on a non-lattice approach, that the the divergence 
is stabilized as the value of N is increased. This is indeed expected since there is a 
meta-stable thermal equilibrium in the theory with a decay rate that is very small 
at large N. 

The numerical results of the internal energy of the thermal gauge theory are 
found in Refs. |32|[33l[3ll[9l[35l[2a[23l[l6l[27]. 

Recently, the authors of Ref. [27] computed the internal energy of the gauge 
theory from Monte Carlo simulations at several lattice spacings and extrapolated 
the result to the continuum limit. Their simulations show that the Polyakov loop is 
largely independent on the lattice spacing for all temperatures probed. The extend 
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of space also experienced very weak lattice effects. However, the internal energy is 
affected by lattice effects. Since the effects of lattice spacing died out linearly in the 
simulations, they were able to extrapolate the energy to zero lattice spacing. They 
used the a' corrected expression Eq. (I13|) while comparing the lattice results with 
that of the gravity theory and found good agreement. 

They also observed that the model became unstable while simulating the theory 
with smaller V, which is related to Hawking radiation in the dual gravitational the- 
oryHHm] Thus one needs larger matrices for a clear comparison with the AdS/CET 
predictions. The results they found were in excellent overall agreement with the 
studies performed in Refs. [32] and [16] . 

In Refs. [THlill] Kadoh and Kamata simulated the topologically twisted version 
of the BESS model. Since they have one additional unintegrated auxiliary field in 
the theory they used a modified expression for the internal energy 

(^{Sb) - - I)) , (96) 

where k denotes the number of unintegrated auxiliary fields, and they simulated the 
theory with k = 1. The Q-exactness of the lattice action implies that the internal 
energy vanishes in the zero temperature limit where the effect of the supersymmetry 
breaking boundary conditions vanishes. 

Their results for the internal energy of the gauge theory coincide with the results 
of the high temperature expansiorP^ at higher values of temperatures (See Eig. 3 
of Ref. [H].) The lattice data smoothly approach the theoretical prediction of the 
gravity side as the temperature is decreased, indicating a strong evidence for the 
gauge-gravity duality. 

Kadoh and Kamata were able to probe the internal energy at the lower tem¬ 
perature region (See Eig. 4 of Ref. [H]). Their investigations show that the lattice 
data approach the gravity prediction and predict that they are likely to coincide 
with it as the temperature decreases further. The temperatures they used in the 
simulations were not low enough to explain the leading behavior on the gravity side, 
unfortunately. To obtain quantitative results for the leading-order term, simulations 
at further low temperatures are required. 

Instead, they studied the contribution of the next-to-leading order term by fitting 
the lattice data using the following formula 

E{T) = 7.41 ® -h C TP, (97) 

where C and p are the fit parameters. From Eq. (US]) , if the duality conjecture is true, 
the value p obtained from lattice data should be 4.6 within the statistical errors. The 
fit performed in Ref. [16] using 5 lattice points within the range 0.375 <T< 0.475 
gave 

C = 9.0(2.6), p = 4.74(35). (98) 

The value of p is consistent with the theoretical prediction from the gravity side 
within about seven percent statistical error. This is, remarkably, the first lattice 
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result of the NLO term, which quantitatively shows the validity of the gauge-gravity 
duality in this thermal system. 

In Ref. [9], the NLO term was estimated from the numerical simulation based 
on the momentum sharp cut-off method using the same fit formula Eq. (EZl) , in a 
little higher temperature region 0.5 < T < 0.7. They obtain the following values for 
C and p after the fit 

0 = 5.55(7), p = 4.58(3). (99) 

They are consistent with the results of Kadoh and Kamata, Eq. (IMl) . within two 
sigma. 

The results in Ref. m were obtained at a finite lattice spacing. It would be in¬ 
teresting to see if the results get better once simulations were performed at different 
lattice spacings and the results are extrapolated to the continuum limit. 

The maximally supersymmetric Yang-Mills quantum mechanics describes the 
open string degrees of freedom attached to the DO-branes, which are decoupled 
from the degrees of freedom in the bulk in the decoupling limit given in Eq. O- 
The above results imply that we can understand the microscopic origin of the black 
hole thermodynamics, including the a' corrections, in terms of the open strings 
attached to the DO-branes. The statement is valid for any value of A and N, and 
lattice is a well suited tool to probe the theory at such arbitrary values of the 
parameters. 


7. The BMN Matrix Model 

Recently, a variation of the BESS matrix model, which has been conjectured to 
describe a discrete light-cone quantization of M-theory on a pp-wave background 
has been formulated.^ This model, known as the BMN model or the Plane Wave 
Matrix Model, can be constructed as a one-parameter deformation of the BESS 
matrix model. 

The dual gravitational geometries to the vacua of the BMN model were con¬ 
structed in Refs. [381139| . These supersymmetric vacuum geometries, including the 
one dual to the trivial vacuum, asymptote to the plane wave solution of M-theory. 

The BMN model has several advantages over the BESS model. It has a discrete 
energy spectrum and a well-defined canonical ensemble. It is convenient to define 
two dimensionless parameters in this model - a dimensionless coupling constant 

5 = 4 - ( 100 ) 

with A the’t Hooft parameter and ft the mass deformation parameter; and 

t=-, (101) 

with T the temperature of the system. These two dimensionless quantities can be 
used to parameterize a two-dimensional phase diagram of the theory. This means 
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that we can use the dual gravitational description at large N and strong coupling 
5 ^ 1 to predict several observables as functions of the dimensionless temperature t. 
The BMN model is expected to have a phase transition whose critical temperature 
should be easy to measure using Monte Carlo simulations. 

In BMN model the flat directions associated with the scalar eigenvalues are re¬ 
moved by the mass terms. Recent efforts to study the BFSS matrix model at strong 
coupling, using both lattice and non-lattice techniques, suffer from the problem of 
flat directions. The existence of these flat directions means that, at finite tempera¬ 
ture, the partition function is formally divergent. It was shown in Ref. [23] that when 
Monte Carlo simulations of this theory are performed, this divergence eventually 
causes the simulations to break down. 

We write down the action of the BMN model 

N / 

^ ^ f Tr ( - Y^iDrXif - + ^'^^rDr'S 

i i<j 

i i—1 i—4 

—2V2iJ€ijkXiXjXk + , ( 102 ) 

where Xi are nine scalars and represents the fermions. The fields of the theory 
are in the adjoint representation of SU{N). Dt is the gauge covariant derivative 
and 

7123 = (103) 

The first four terms of Eq. (11021) are the same as that of the BESS matrix model. 
The rest of the terms give supersymmetry preserving masses to the fields and add 
a Myers term. As mentioned above, the addition of these terms gives this model 
technical advantages compared with the BESS. Firstly, in the limit /i —>■ oo the 
model becomes weakly coupled and can be studied using perturbative methods .1^ 
The addition of these terms lifts the moduli space of the BFSS model giving instead 
a discrete set of vacua.^ This is indeed important for numerical simulations as it 
means the theory is well defined in the canonical ensemble. 

If we consider the scalar fields Ai for z = 1, 2, 3 we can write the potential as 


V = -Tr 


fiXi 




k \Xi , Xj 


(104) 


which is minimized when the scalar fields Xi are proportional to generators of 
SU (2). The discrete vacua correspond to the various ways of forming NxN matrices 
that generate SU(2), and can be put in one-to-one correspondence with the integer 
partitions of TV. It was shown in Ref. [41] that the vacua are free from quantum 
corrections and thus persist at strong coupling. 
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To write down the lattice action, the parameters of the continuum BMN action 
Eq. (II02I) are rescaled in the following way 


V ^ nMX, 

(105) 

Dt — y iJiA4^ 

(106) 


(107) 

4' (^M)^4', 

(108) 


where M is a dimensionless number that will become the number of lattice sites. 
The lattice action is 

5 ^ f dr Tr ( - 

i i<j 

+ Y, >P| - ^(X,)^ 

i i—1 i—4 

-2V2^e.,kX,X,X, + AvI/^7,23'1') • (109) 

All the fields in the above action are dimensionless. The next step is to integrate 
out the fermions to get the following partition function 

Z = J DADX Pf{M)e-^^, ( 110 ) 


with the bosonic action 


Sb = 




2A 


f 


dr Tr - 




i<j 


4M2 


^(V,)2 - 2V2-^e,jkX,XjXk), 


i—4 


1 

AP 


3 




(111) 


and the fermion operator 


M =-frDr+J2-/i[Xi, •] + ^ 7123 - 


The lattice theory has an overall coupling 

/V/3 


K = 


2A 




( 112 ) 


(113) 


where (3 is the inverse temperature. 

One could then proceed to simulate this model using RHMC algorithm. The 
continuum limit of the theory is defined by taking the number of lattice points 
M —>• 00 with ^^V/3/(2A) fixed. 

In Ref. [32] Catterall and van Anders presented the first results of the lattice 
simulations of the BMN model. They used a naive lattice action, and concentrated 
on studying the Hagedorn/deconfinement transition in the model. They studied 
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the theory at a fixed temperature, measured in units of the mass deformation /r, 
as a function of the ’t Hooft conpling A (measured in the same nnits). They used 
a quenched approximation to explore the dependence of the critical behavior on 
the rank of the gauge group N and the number of lattice points M. The quenched 
approximation is much less computationally demanding and can be used to estimate 
reasonable choices of parameters for simulating the full theory. Their simulations 
indicate that one can get a reasonable approximation of the continuum large N 
behavior with modest values of both N and size of the lattice. By simulating the 
model at fixed temperature over a range of coupling they showed that the model 
exhibits a deconfinement transition when the’t Hooft coupling is of order one. They 
also found that the critical value of the coupling in the case of quenched simulations 
was of the same order as in the theory with dynamical fermions. 

According to the conventions of Ref. [43] , Catterall and van Anders simulated 
the BMN at fixed finite temperature, T/p. = 1/3, and as they varied the coupling, 
they observed a first order phase transition for 0.03 < 5 < 0.045. This result is not 
in direct contradiction with the results of Ref. [43] but it implies a non-monotonic 
behavior of the critical temperature as a function of the coupling g, complicating 
the phase diagram computed in Ref. [43]. 

In Ref. [43] the authors found a Hawking-Page like phase transition in the dual 
gravitational description of the BMN model and predicted the strong coupling limit 
of the critical temperature 


lim = 0.105905(57). (114) 

g-J-oo g 


It would be remarkable to confirm this prediction with Monte Carlo simulations of 
the BMN model at strong coupling. 


8. Dl-branes and Black Hole - Black String Transitions 

In this section we briefly discuss the lattice results from gauge-gravity duality con¬ 
necting two-dimensional maximally supersymmetric Yang-Mills theory and type HB 
supergravity solution describing the thermal vacuum of black hole carrying electric 
Ul-brane charge. 

On the gauge theory side we have a large N finite temperature two-dimensional 
maximally supersymmetric SU(N) Yang-Mills theory, in the ’t Hooft limit, with 
coupling A = NgYM-: with the spatial direction compactified. In Euclidean time, this 
implies the Yang-Mills theory is defined on a rectangular 2-torus, with time cycle 
size /3, and space cycle size R. The fermion boundary conditions distinguish the two 
cycles, being anti-periodic on the time cycle and periodic boundary conditions on 
the space cycle. 
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The action of the Euclidean theory is 

5'= y J ^drdx Tv + 

+ fermions^, (115) 

where i, j = 1,..., 8 and (j)i are the 8 adjoint scalars, and t is the coordinate on the 
time circle, and x the coordinate on the space circle. 

It is convenient introduce two dimensionless couplings in the theory, 

Tr = A^/3, rx=\^R, (116) 


which give the dimensionless radii of the time and space circles respectively. 

We will be interested in the expectation values of the Polyakov loops on the time 
and space circles. 


{Pt,x) 


1 

V 


Tr 


exp 



(117) 


At large V, they give order parameters for confinement/deconfinement (or center 
symmetry breaking) phase transitions in the thermal gauge theory. 

There are several interesting limits of this theory, as discussed in Refs. IHISI]. 
We are interested in the large torus limit, I <C rx,rr where the string theory dual 
may be described by a supergravity theory.E^ Having a supergravity description 
of the full string theory dual allows certain behaviors of the theory to be studied 
using simple semi-classical gravity reasoning, which allows powerful predictions to 
be inferred for the dual SYM. 

The dual type IIB string theory is given by the ‘decoupling limit’ of N coincident 
ZH-branes.l^ Since our Euclidean SYM is defined on a torus, the string dual is too, 
being at finite temperature and having one spatial direction compactified into a 
circle radius R with periodic fermion boundary conditions. 

One finds that for 1 <C r,- <C this string theory can be described effectively by 
its supergravity sector. The type IIB supergravity solution describing the thermal 
vacuum is a black hole, carrying electric Z?l-brane charge. In the type IIB regime, 
1 Tt- <C we expect Pr ^ 0 but = 0. In the type BA regime, where I <C r,- 
and rr, we have Pt ^ 0, and Px ^ ^ for < Ccrit?'r and Px = 0 for 

> Ccritfr, with Ccrit an order one constant with Ccut > 2.29. We note that in 
the regime where both type BA and IIB apply, they give consistent results. Thus 
in the large torus, supergravity regimes, the SYM is always deconfined in the time 
direction, and there is a first order deconfinement/confinement transition in the 
space direction at 


X — Ccrit^r ■ 


(118) 


In the simplest picture, the above relation characterizes the Gregory-Laflamme (GL) 
type first order phase transitiorP^ for I <C r^, with Ccut > 2.29. 
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8.1. Supersymmetric Lattices for 2D SYM 


In this section we briefly discuss the manifestly supersymmetric lattice action used 
to simulate the strongly coupled two-dimensional SYM theory. The lattice theory 
is based on a twisted version of the four-dimensional maximally supersymmetric 
theory. The continuum twist of A/" = 4 we are interested in was first written down by 
Marcus.l^ As shown by Kaplan and UnsaP^^this twisted version of the theory can 
be recovered by taking the continuum limit of a lattice theory constructed using the 
idea of orbifold projection. CatteralP^ has constructed a direct lattice discretization 
of this theory combining the ideas of topological twisting and Kahler-Dirac fermions. 
These two lattice formulations lead to identical supersymmetric lattices for this 
theory.This four-dimensional twisted theory can be most compactly expressed as 
dimensional reduction of a five-dimensional theory In which the ten bosonic fields 
(one four-component gauge field and six scalars) are realized as the components of a 
complexified five-dimensional gauge field Am, m = 1 ... 5, while the sixteen single¬ 
component twisted fermions are realized as the sixteen components of a Kahler- 
Dirac field {rj, ipm, Xmn}- Twisting leads to a scalar fermion in the theory - it implies 
the existence of a nilpotent supersymmetry which will be preserved in the lattice 
theory. 

The scalar supersymmetry charge acts on the continuum fields In the following 
way 


QAm — 

(119) 

QV'm = 0 

(120) 

QAm=0 

(121) 

QX^n — Amn 

(122) 

Qt] = d 

(123) 

o 

II 

0 ^ 

(124) 


The field Am is the complex conjugate of Am and J-mn is the complexified curvature 
made out of Am- The scalar field d Is an auxiliary field that Is Included to close 
the Q supersymmetry algebra and Is subsequently Integrated out of the final lattice 
action. 

The action of the twisted theory can be written as the sum of two terms - a 
Q-exact piece and a Q-closed term. 

The Q-exact piece has the form 


S = 


9ym 


Q I Tr ( Xm-nA"mn 4“ y[Dm,‘Dm] 


and the Q-closed term 


•^closed — 


1 


Tr emnpqrXgrDpXn 


(125) 


(126) 


^9ym J 

The Q-exact piece is trivially Q supersymmetry invariant. The supersymmetric 
invariance of the Q-closed term can be shown using the Bianchi identity for the 
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complexified covariant derivatives 

^mnpqr'^p^ qr — 0 - 

We can discretize the theory in a straightforward manner. The complex contin¬ 
uum gauge fields are represented as complexified Wilson gauge links 

Uf,in) = (128) 

living on positively oriented links e^, /r = 1,..., 4, of a four-dimensional hypercubic 
lattice with integer sites n. The field is placed on the body diagonal of the 
hypercube corresponding to a relative position vector 

65 = (-1,-1,-1,-1). (129) 

Notice that we have the restriction 

5 

^ = 0 (130) 

m—1 

in order for the action to be gauge-invariant on the lattice. These fields transform 
in the usual way under the U{N) lattice gauge transformations. Eor example, we 
have the following gauge transformation rule for the complexified gauge link 

Um{n) ^ G{n)Um{n)G'^{n-\-em), (131) 

for G e UiN). 

Supersymmetric invariance then implies that '0m(n) live on the corresponding 
links and transform identically to Wm(n). The scalar fermion r]{n) is associated with 
a site and transforms like a site field 

r]{n) G'(n)? 7 (n)G^(n). (132) 

The field Xmn is placed on the diagonal links connecting sites n-|-em + e'„ and n. 
This choice of orientation will again be necessary to ensure gauge-invariance. The 
scalar lattice supersymmetry transformation is identical to that in the continuum 
after the replacement 

(133) 

Most importantly it remains nilpotent, which means that we can guarantee invari¬ 
ance of the Q-exact part of the lattice action by replacing the continuum fields by 
their lattice counterparts. 

We also need a prescription for replacing continuum derivative operators 
by gauge-covariant finite difference operators. We use the following prescrip¬ 
tions A forward lattice covariant difference operator is used to denote a 

curl-like operation on the twisted fields 


T>^Vn(n) =Wm(n)/„(n-Fem) - /„(n)Wm(n -h e„). 


(134) 
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while backward lattice covariant difference operator is used to denote a divergence¬ 
like operation 

Vm(n) = /m(n)Wm(n) - Wm(n - em)/m(n - e^). (135) 


They reduce to the usual adjoint covariant derivatives in the naive continuum 
limit. It is also guaranteed that the resultant expressions transform covariantly 
under lattice gauge transformations. 

It is possible to discretize Eq. (11261) in such a way that it is indeed exactly 
invariant under the twisted supersymmetry 


»5*closed — 


Sg 


YM 


■ETr 


^mnpqr\.qr{.^ + 671 + Cp)'^p Xmn(n + 6p), ( 136 ) 


and can be seen to be supersymmetric since the lattice field strength satisfies an 
exact Bianchi identity 


^mnpqr'^p ^qr — 0 


(137) 


Assembling all these together we arrive at the supersymmetric lattice action 
(See Refs. [49l[52] for reviews) 

1 . . 1 f^(-) ^2 


5 = ^closed + ^ ETr (-FL(n)^mn(n) -b - (v^ 

9 ym ^ ^ ^ ^ 


-X™„(n)T'[+V„](n) - ?7(n)T>)„ Vm(n)), 


■(-) 


(138) 


where we have taken the Q-variation and the auxiliary field has been integrated 
out. This action is gauge-invariant, free of doublers and possesses the one exact 
supersymmetry charge Q. 

In order to arrive at the two-dimensional theory we perform a simple dimen¬ 
sional reduction along two lattice directions using periodic boundary conditions. 
The resultant lattice action corresponds in the naive continuum limit to the target 
sixteen supercharge Yang-Mills theory in two dimensions. In this limit its exact 
supersymmetry is enhanced to four supercharges corresponding to the four scalar 
fermions that now appear in the dimensionally reduced theory.l^ 

We will be interested in this theory at the large N limit with’t Hooft coupling 
A. The lattice theory is then governed by the coupling 


K = 


NLM 

2r2 


(139) 


where L and M denote the number of lattice sites in the spatial and temporal 
directions, respectively. 

We use periodic boundary conditions for the fields on the remaining spatial 
circle and anti-periodic boundary conditions for fermions in the temporal direction 
in order to access the theory at finite temperature. Simulations can be carried out 
using RHMC algorithm. It has been shown that the existence of a non-compact 
moduli space in the theory may render the thermal partition function divergent.!^ 
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In order to regulate this divergence one introduces mass terms for the scalar fields 
appearing in the lattice action with a dimensionless mass parameter 

m = TOphys/3- (140) 

The mass terms added to the action have the form 

2 

= (ulU^ + {Up^y^ - 2 ) . (I4I) 

9ym „ ^ ^ 

Such mass terms are effective at suppressing arbitrarily large fluctuations of the 
exponentiated scalar fields and reduces to simple mass terms for small fluctuations 
characterizing the continuum limit. These infrared regulator terms break super- 
symmetry softly and lift the quantum moduli space of the theory. In Ref. [53l [54] 
simulations were performed for a range of the mass parameter in order to allow for 
an extrapolation m —>■ 0. 


8 . 2 . Results on Gregory-Laflamme Phase Transition 


In this section we present the numerical results of the studies performed in Refs. 
[551151] . The simulations are based on the Polyakov lines for both the thermal and 
spatial circle. These are defined in the usual way 


Px 

Pr 


(iTrnyyu^j), 

(142) 

(pni-yuap, 

(143) 


where the unitary piece of the complexified link Up, has been extracted to com¬ 
pute these expressions. The spatial and temporal Polyakov lines are evaluated as 
functions of r^. for various lattice sizes and number of colors. 

The temporal Polyakov line remained close to unity in the simulations over a 
wide range of r^, indicating that the theory is (temporally) deconfined. The spatial 
Polyakov line has a different behavior, they take values close to unity for small Tt 
and fall rapidly to plateau at much smaller values for large r,-. It is tempting to 
see the rather rapid crossover around Tt ^ 0.2 as a signal for a would be thermal 
phase transition as the number of colors is increased. This conjecture is seen to be 
consistent with the lattice data. The plateau evident at large falls with increasing 
N and the crossover sharpens. This is consistent with the thermal system developing 
a sharp phase transition in the large N limit. 

Simulations show that the m = 0 model does exhibit the same thermal instability 
observed in the case of sixteen supercharge quantum mechanics for sufficiently low 
temperature r^- >> I, in agreement with the general arguments given in Ref. [23] . 

The gravity prediction for the parametric behavior = Ccrit?'r is consistent with 
the lattice data, and Ccrit was estimated to b^ 


Ccrit — 3.5, 


(144) 


which indeed obeys the gravity prediction that Ccrit is order one and Ccrit > 2.29. 
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The value of the ratio 


a = 


Ccrit 

2 ^ 


( 145 ) 


gives the ratio of the GL thermal phase transition temperature to the GL dynamical 
instability temperature (the minimum temperature to which uniform strings can be 
supercooled), so 


a = 


^GL phase 
Tgl instab 


(146) 


The GL instability temperature is knowrP^ (corresponding to the behavior = 
2.29rT- at strong coupling) however, the GL phase transition temperature is not 
known in the gravity theory as the localized solutions have not been constructed 
yet. 

The lattice estimatiorP^ 


a-1.5 (147) 

provides a prediction for the thermal behavior of the gravity solutions. This was the 
first time a prediction about the properties of non-trivial classical gravity solutions 
had been made from the Yang-Mills side of a gauge-gravity correspondence. Since 
the dual localized black hole solutions have not been constructed, this constitutes a 
prediction for these non-trivial gravity solutions, which hopefully will be tested by 
their construction in the near future. 


9. Future Directions 

In this review we have focused on maximally supersymmetric lattice Yang-Mills 
theories in one and two dimensions and showed how they can validate and provide 
insight into the gauge-gravity duality conjecture. For the case of one-dimensional 
theories we have detailed the duality conjecture connecting DO-branes and black 
hole geometries at finite temperature to maximally supersymmetric thermal gauge 
theories at the same temperature. The internal energy of a quantum black hole 
with its leading order quantum gravity and string theory corrections has been de¬ 
rived recently in the gravitational theory. The predictions from supergravity can 
be verified by simulating the dual strongly coupled thermal gauge theories on the 
lattice and calculating the internal energy of the thermal system for a given number 
of colors and’t Hooft coupling. Recent simulations indicate that the gauge theory 
results give excellent agreement with the quantum gravity predictions, thus proving 
strong evidence for the validity of the duality conjecture. We detailed the naive 
and manifestly supersymmetric constructions of maximally supersymmetric lattice 
Yang-Mills theories and the simulation algorithms used to probe the strong cou¬ 
pling regimes of these gauge theories at finite temperature and number of colors. 
We have also briefly detailed the results on the duality in two-dimensional sixteen 
supercharge thermal gauge theories and corresponding supergravity solutions. We 
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also note that the existing results for the two-dimensional system can be improved 
using improved lattice simulation methods recently proposed in Refs. |551156j . 

It would be interesting to extend the analysis of the gauge-gravity duality to 
three- and four-dimensional Yang-Mills systems, which are thought to be dual to 
D2- and D3-branes systems. Such theories are easily accessible through the exact 
lattice supersymmetry formulation. 

Another interesting matrix model to investigate in the context of gauge-gravity 
duality is the Berkooz-Douglas matrix model.^ This model is holographically dual 
to the back-reacted D0-D4 system. Investigations on this model can give more ex¬ 
tensive tests of the gauge-gravity duality. 

Numerical simulations of the sixteen supercharge quantum mechanics model 
become increasingly difficult for smaller N and lower temperatures. The simulations 
suffer from technical difficulties related to the divergence of the eigenvalues of the 
bosonic matrices while sampling the important field configurations that contribute 
to the action. It would be remarkable if we could tackle this problem using new 
numerical or analytical methods to provide more insight into the gauge-gravity 
duality and the internal structure of quantum black holes. 
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